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ABSTRACT 

Triaxial neutron stars may be important sources of gravitational radiation for 
the forthcoming generation of interferometric gravitational wave detectors such as 
LIGO, VIRGO, and GEO600. We investigate the viscosity triggered bar mode secular 
instability of rapidly rotating neutron stars by means of a perturbation analysis of 
numerically constructed "exact" general relativistic axisymmetric star models. In 
the theoretical approach, only the dominant parts of the nonaxisymmetric terms 
of the 3D-Einstein equations are taken into account. A comparison of our results 
with previous studies of Newtonian polytropic stars confirms James' classical result 
7crit = 2.238 for the critical polytropic index. Beyond the Newtonian regime, 7crit 
reveals a slight increase toward highly relativistic configurations. Six out of twelve 
employed realistic dense matter equations of state admit the spontaneous symmetry 
breaking for masses above 1.6 M©. 



1. Introduction 

Rapidly rotating neutron stars are highly relativistic objects, and, provided there is some 
physical process operative, which induces a significant deviation from axisymmetry, may be 
important sources of gravitational radiation. A transition toward triaxial configurations can 




Fig. 1. — Triaxial maximum rotation 
Newtonian polytropic star for 7 = 2.4. The 
ellipticity in the equatorial plane is about 
0.1. Notice the cusp at the stellar equa- 
tor in direction of the semimajor axis {x- 
axis), where the rotation is Keplerian, and 
which is absent along the semiminor axis 
(y-axis). 



occur, when the ratio T/|iy| of rotational kinetic energy and gravitational potential energy 
reaches some critical value^'^-*. The equation of state, hereafter EOS, has an important influ- 
ence on the development of a triaxial instability since the neutron star matter must be stiff 
enough to admit a maximum angular velocity Qk higher than the critical one. Homogeneous, 
incompressible fluid bodies, rotating at moderate constant angular velocity, take the shape of 
some oblate axisymmetric Maclaurin spheroid. Maclaurin spheroids are dynamically unstable 
for T/|iy| > 0.2738, but there are two other families of triaxial ellipsoids which bifurcate from 
the Maclaurin sequence earlier at T/|iy| =0.1375. We focus on the sequence of Jacobi ellipsoids 
which rotate uniformly about their smallest axis in an inertial frame. The evolution toward 
a Jacobi ellipsoid at the bifurcation point is triggered by some viscous dissipative mechanism. 
This instability is called a secular one, because it evolves on the associated viscous timescale 
which is much longer than the dynamical one. The astrophysical realization corresponds to 
the scenario of a cold highly viscous binary neutron star, being spun up by accretion from its 
companion. Viscosity dissipates mechanical energy, while it preserves angular momentum. As 
a consequence, the Maclaurin spheroid, once the instabihty point is reached, develops toward 
a Jacobi ellipsoid which exhibits the lowest rotational kinetic energy for a fixed value of the 
angular momentum. In this final state, viscous dissipation has ceased, and the star is rotating 
rigidly. Former studies had been performed at the Newtonian^'^'^'^-* or Post-Newtonian''''^'' level 
which is by no means adequate for highly relativistic objects like neutron stars. Furthermore, 
the neutron star matter was modeled by a simplifying polytropic equation of state. 



2. Numerical models of triaxial neutron stars 

Before the symmetry breaking sets in, the neutron star, modeled as a rigidly rotating per- 
fect fluid, can be considered as stationary and axisymmetric. Making the further assumption 
that no meridional matter currents are present, a favourable choice of elliptic field equations 
was given by Bonazzola et al.^^ The numerical code^-*, based upon these equations, relies on a 
spectral method, and allows to compute neutron star models with a precision of ~ 10~^^ in the 
spherical symmetric case and ~10~^ for maximum rotation configurations, when a 7 = 2 poly- 
tropic, analytic EOS is used, as well as ~10~^ for realistic EOS^°^ The neutron star models are 
"exact" in the sense that the full Einstein equations are solved without any analytic approxima- 
tion, while the numerical integration covers all space and respects the exact fiat space boundary 
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Fig. 2. — Critical polytropic index 
7crit as function of the lapse func- 
tion Nc, measured at the centre of 
the star. Black boxes indicate con- 
figurations unstable with respect to 
radial oscillations. 
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condition at spatial infinity. When the symmetry breaking occurs, spacetime is neither sta- 
tionary nor axisymmetric. However, at the very beginning, the deviation from axisymmetry 
is sufficiently small, and the emitted gravitational radiation may be neglected. Under the ad- 
ditional assumption of rigid rotation, the helical symmetry of spacetime is preserved. Taking 
into account only the dominant non-axisymmetric terms, one basically recovers the original 
equations of the axisymmetric case^^-*. The only difference is that the lapse function — its 
logarithm u = log N reduces to the Newtonian potential in the weak field approximation — 
and the matter fields are three-dimensional quantities which depend on r, 6, and ip = (p — Qt, 
where fl is the angular velocity of the star. As a consequence, the present approach is exact 
for the stationary and axisymmetric relativistic, as well as for the fully three-dimensional New- 
tonian case. After relaxation to a particular axisymmetric configuration, a small perturbation 
= eHc{r sin 6 cos tp)"^, which excites the 1 = 2, m = ±2 bar mode, is added to v = \ogN, 
and the growth of 5v during the subsequent relaxation is followed. He denotes the central 
log-enthalpy, and e is a small parameter of the order of 10^^. A particular configuration is 
secularly unstable if the perturbation increases, and the fiuid body evolves subsequently toward 
a triaxial spheroid. It is conditionally secularly stable if the perturbation tends to zero. In this 
case, however, three-dimensional terms of higher relativistic order, which have not been taken 
into account at the present level of approximation, may induce the symmetry breaking. 



3. Results for polytropic and realistic equations of state 

The investigation of the symmetry breaking of relativistic polytropic stars avoids the prob- 
lems associated with the use of realistic dense matter EOS^^-', and allowed us to perform some 
comparison with former results obtained for Newtonian stars'*'^'^-'. In particular, we have con- 
firmed James'^-* value 7crit = 2.238 of the critical polytropic index, where secular instability sets 
in at maximum angular velocity, as well as the asymptotic ratio T/|iy|^^oo = 0.1375. No results 
in the fully relativistic regime had been reported so far. In the Newtonian case, polytropic stars 
obey a scaling law, and the critical polytropic index is a global constant. Relativistic effects, 
however, are supposed to influence the symmetry breaking. Fig. |^ shows the dependence of 7crit 
on the central value of the lapse function Nc which is an appropriate signature of the relativistic 
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Tab. 1. — Neutron star prop- 
erties according to various 
EOS: M^^* is the maximum 
mass for static configurations, 
-Mmax is the maximum mass 
for rotating stationary con- 
figurations, Pk is the cor- 
responding Kepler ian period, 
Phreak is the rotation pe- 
riod below which the symme- 
try breaking occurs, i^c,brcak 
is the central log-enthalpy 
at the bifurcation point and 
Mbreak is the Corresponding 
gravitational mass. The EOS 
are ordered by decreasing val- 
ues of M±^l 



character of the star. After a very slight decrease in the weakly relativistic regime it grows by 
about 10 % for strongly relativistic configurations. One concludes furthermore that any neutron 
star, built upon a polytropic EOS and stable against radial perturbations, becomes secularly 
unstable for 7 > 2.33 at some critical angular velocity ficrit<^K- The results for twelve realistic 
EOS of neutron star matter are presented in Tab. |I]. A description of each EOS can be found 
in Salgado et al.^°\ where they had already been employed to construct realistic high precision 
neutron star models. Six equations turn out to be stable up to the mass-shedding limit ^k- 
Only the stiffest ones allow for the symmetry breaking. The correlation with the corresponding 
maximum mass is not strict because of the density dependence of the EOS. This is notably 
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Fig. 3. — Measured masses 
of 18 neutron stars. Ob- 
jects in massive X-ray bina- 
ries are at the top, radio pul- 
sars and their companions at 
the bottom. Data compiled 
by Thorsett et al.^^^ is marked 
by filled circles. The empty 
circle indicates the mass of 
the recently discovered ra- 
dio pulsar PSR J1012-h5307, 
measured by van Kerkwijk et 
Si\}^\ The continuous line 
indicates the estimate of the 
mean value of the neutron star 
mass distribution derived by 
Thorsett et al. and the dashed 
lines the corresponding la in- 
terval. 
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the case for the WGW EOS which is the stiffest one at moderate densities. For the present 
sample of reahstic EOS, Mbreak has a lower bound of 1.645 Mq. This is a low value compared 
with maximum masses of rotating neutron stars built upon a stiff EOS, e.g. in Tab. |I|. Un- 
fortunately, the real distribution of neutron star masses is still quite uncertain. Fig. ^ shows 
the measured masses of 17 neutron stars, compiled by Thorsett et al}^\ and yielding an es- 
timate of the average neutron star mass of Mns = 1.35 ± 0.27 M©, as well as a recent mass 
measurement of PSR J1012+5307 of 2.35 ± 0.434 Mq^^). At least the masses of the two X-ray 
binary neutron stars 4U 1700-37, Vela X-1, and that of the radio pulsar PSR J1012+5307 are 
well above our minimum mass of 1.645 M©. If they prove to be valid after further reduction 
of the observational errors, we may expect a significant number of neutron stars susceptible of 
encountering the spontaneous symmetry breaking. Their observation would provide us with 
useful information about the generation of gravitational waves and constraints on real neutron 
star matter EOS. 
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ETOILES A NEUTRONS TRIAXIALES — UNE SOURCE 
EVENTUELLE DE RAYONNEMENT GRAVITATIONNEL 

Des etoiles a neutrons triaxiales pourraient etre des sources de rayonnement gravitationnel 
importantes pour la prochaine generation de detecteurs interferometriques d'ondes gravitation- 
nelles comme LIGO, VIRGO ou GEO600. Nous etudions I'instabilite seculaire du mode barre des 
etoiles a neutrons en rotation rapide par I'intermediaire d'une analyse perturbative de modeles 
stellaires axisymetriques "exacts" en relativite generale. Dans I'approche theorique on ne tient 
compte que de la partie dominante des termes non-axisymetriques des equations d'Einstein. 
Une comparaison avec des etudes precedentes d'etoiles newtoniennes confirme le resultat clas- 
sique de James a savoir que I'indice polytropique critique vaut 7crit = 2.238. Au-dela du regime 



newtonien 7crit montre une faible croissance vers les configurations fortement relativistes. S 
sur douze equations d'etat de la matiere dense utilisees admettent la brisure spontanee de 
symetrie pour des masses superieures a 1.6 Mq. 



